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Abstract 

We discuss how one can determine the average kinetic energy of the heavy quark 
inside heavy mesons from differential distributions in B decays. A new, so-called 
third, sum rule for the h — > c transition is derived in the small velocity (SV) limit. 
Using this sum rule and the existing data on the momentum dependence in the 
B — ► D* transition (the slope of the Isgur-Wise function) we obtain a new lower 
bound on the parameter /j% = (2Mb)" 1 (B\b(iD) 2 b\B) proportional to the average 
kinetic energy of b quark inside B meson. The existing data suggest ft 2 > 0.4 GeV 2 
and (from the "optical" sum rule) A > 500 MeV, albeit with some numerical uncer- 
tainties. 
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1. In two recent papers jIL 0] it was shown how the operator product expansion 
(OPE) allows one to derive various useful sum rules for heavy flavor transitions in 
the small velocity (SV) limit ||- Non-perturbative corrections are included in the 
theoretical side of the sum rules in the form of an expansion in the inverse powers 
of the heavy quark mass. In Ref. 0] the so called first sum rule at zero recoil was 
obtained which was then used for estimating the deviation of the B —* D* transition 
form factor from unity at zero recoil to order 0(Aq CD ). Another sum rule analyzed 
in |L| yields a field-theoretic proof of the inequality 

/4 > /4 (i) 

where /i^ and p? G are related to the kinetic energy and the chromomagnetic operators, 

A = ^- B {B\l{Wfb\B), nl = ^L-(B\b(i/2)*Gb\B). (2) 

(Previously this inequality was obtained within a quantum-mechanical approach 
J|, HO In the present paper we exploit similar ideas to get a new sum rule in 
the SV limit which relates fx\ to the average square of the excitation energy of the 
final hadronic state X c in the B — > X c semileptonic transitions. At present the 
corresponding inclusive differential distribution is not yet measured. However, we 
use the existing data on5-> D*lv decays near zero recoil to get a lower bound on 
/i^, with no reference to \j? G . The bound involves the slope of the Isgur-Wise function 
extracted from the momentum dependence of the B — > D* transition. Numerically 
this bound turns out to be close to that of Eq. ([[]). 

2. The general method allowing one to derive the sum rules in the SV limit 
is presented in Ref. [|l[. Here we remind only some basic points primarily for the 
purpose of introducing relevant notations. Operator product expansion is applied 
to the transition operator || [7j 

f ab (q)=tJd i xe^T{jl(x) 3b (0)} (3) 

where j a denotes a current of the type cT a b with an arbitrary Dirac matrix r a ; q is 
the momentum carried away by the lepton pair. The average of T ab over the heavy 
hadron state H b with momentum pn b represents a forward scattering amplitude (the 
so-called hadronic tensor), 

hab(PH b ,q) = ^7l—( H b\Tab\H b ). (4) 

The observable distributions are expressed through the structure functions w ab , 

Wab = (l/i) disc h a b- 



In the HQET limit ||, when one neglects l/m&, l/m c corrections, the hadronic 
tensor h a h is defined by one invariant function h for any matrix T a in the current j a , 
namely: 



hah — Cahh , C a h — Tr 



i a ~ 1 6 



2 2 

Here T a = 701^70 and f 1^, t^ are 4- velocities of the initial and final hadrons, 



(5) 



^ - MhT ' ^ - M Hc ' (6) 

(M# 6 and Mh c can be substituted by m and m c , correspondingly, in the leading 
approximation). The invariant function h depends on two scalar invariants available 
in the process, namely (viq) and q 2 . In what follows we will assume that the hadron 
H D is at rest; the first invariant then reduces to go- Moreover, in studying the 
transitions h — > c at zero recoil or in the small velocity (SV) limit, it is convenient 
to employ directly the spacelike momentum transfer q 2 = (v\q) 2 — q 2 as the second 
argument of h. 

Taking into account higher order in l/mf, iC we loose, generally speaking, this 
property of the factorization of h a b into a universal kinematical structure times one 
hadronic function h. In the general case, the hadronic tensor can be decomposed in 
terms of possible covariants [0 (their number depends on the Lorentz structure of 
the currents) with coefficients hi. In particular, in the case of the vector and axial- 
vector currents we deal with the functions hf v , hf A and hY A , i = 1, ..., 5 introduced 
in Ref. 0. They are independent functions. However, in the leading order they all 
are expressible in terms of h, namely, 



h 



h_ - ^L^hAA _ h AA. 

1 + v iv 2 2 m b 



h\ *- m cuVV uVV 

-ho = m r h 



and 



-2 — "Wb ) 

v x v 2 2 m h 



h = —m c h^ . 



The functions hi not listed here vanish in this approximation. The expressions for 
all functions hi up to order 1/m 2 can be found in 0. 

Although the universal factorization above is not valid for all non-perturbative 
corrections it still holds for those corrections that are relevant for the third sum rule 
to be derived below. We will explain this point shortly. Since it is not important 
what hadronic function we deal with - they all lead to one and the same third sum 
rule - we will use h AA in our derivation. Thus, we consider the transitions of the B 
meson induced by the axial- vector current, 

A" = C7/.75& • 



To single out h AA one must consider the spatial components of the axial current 
generating the transitions of the B meson to D* and the corresponding excitations. 

In the sum rules at zero recoil (q = 0) were obtained; here we will work at 
small but non-vanishing values of \q\. The terms 0{(f) will be kept while those of 
higher order in |g| will be neglected. 

To start the derivation we consider h AA (qo,q) in the complex q plane (q is 
assumed to be fixed, and Aq C d "C |<f| <C M d ). More exactly, let us shift q by 
introducing the quantity 

£ — %max — Qo (7) 

where 

q 0max = M B - E D , , E D . = M D * + -f—. (8) 

When e is real and positive we are on the cut. The imaginary part of hi is given 
by the "elastic" contribution of D* plus inelastic excitations. For what follows it is 
crucial that all these contributions are positive definite. 

For negative e we are below the cut, and the amplitude h AA can be computed - 
and it actually was || [H| - as an expansion in Aqcd/^c,^ For our purposes it is 



sufficient to limit ourselves to the correction terms of the first and the second order 



in Aqcd- This is exactly the approximation adopted in || [T(| [ITJ, and expressions 
obtained there will be used below. 

At the next stage we assume that Aq CD <C |e| <C m bjC . The amplitude h AA 
is expanded in powers of Aq CD /e and ejm^^. Polynomials in e can be discarded 
since they have no imaginary part. We are interested only in negative powers of 
e. The coefficients in front of l/e n are related, through dispersion relations, to the 
integrals over the imaginary part of h AA with the weight functions proportional to 
the excitation energy to the power n — 1. Thus, the first sum rule considered in 
Ref. corresponds to n = 1; the second sum rule (sometimes called optical or 



Voloshin's sum rule [12 1 , see also [13, 14|) corresponds to n = 2. The lower bound 



on fi^ - our main aim in this work - stems from the third sum rule, i.e. we need to 
analyze the coefficient in front of 1/e 3 in the expansion of h\. 

The 1/e expansion can be read off from Eq. (A.l) in Ref. 0. One technical 
element of the derivation deserves a comment. The theoretical expression for the 
amplitude h AA presented in knows nothing, of course, about the meson masses; 
it contains only the quark masses. Correspondingly, it is convenient to built first 
the expansion of h AA in an auxiliary quantity, 

q 2 
e q = m b - E c - q , E c = m c + - — . (9) 

2m c 

Then, if necessary, we reexpress the expansion obtained in this way in terms of e. 
The difference between e q and e is C(Aq CD • (f/ml c ) and 0(AQ CB /m bc ). It will 
be seen shortly that to our accuracy this difference can be simply ignored in the 
third sum rule in the SV limit. It can not be discarded, however, in the second sum 



rule. (The situation is quite different from that which took place in the sum rules 
at zero recoil, see 0. At zero recoil the difference between e and e q was absolutely 
important.) 

The expression for h AA in Eq. (A.l) in || has the form 

-h AA = [( mb+mc - qo ) + 0(Al iCD /m b )]-+0(A 2 QCD )- + -(m b +m c -q )fxlf- (10) 

Z Z o Z 

where 

z = e q (2E c + e q ). (11) 

Notice the similarity of the coefficient in front of 1/z 3 and the leading part of the 
coefficient in front of 1/z. This is not accidental. The terms 1/z 3 appear only as the 
expansion of the second order in nq of the denominator (m^v — q) 2 — 2qir (see Ref. 
||) and, therefore, preserve the same universal factorization which was pointed out 
above in the HQET limit. 

Expanding in e q /2E c we observe that 1/ z n reduces to 1/e™ plus all lower powers 
of l/e q plus a polynomial in e q . The next step is eliminating e q in favor of e. The 
term 1/z 3 comes with a coefficient fi 2 ■ q 2 ; hence here the difference between 1/e and 
l/e q is of higher order and can be neglected. By the same token to order C(Aq CD ) 
one can substitute l/e q by 1/e in 1/z 2 . As far as 1/z is concerned here we must 
reexpress l/e q in terms of 1/e, 



1 (r _,- 

2 



+ ^-^ + ... (12) 



e q e e 

Next terms in Eq. dl2|) are irrelevant since they lead to corrections of higher order 
in Aqcd and/or \q\. This observation is crucial since it tells us that the 1/z part 
contributes only to the first and the second sum rules; it generates no 1/e 3 terms. 
As a result h AA has the form 

~h AA = \ [l - £j + 0(A 2 QCD K)) + i (0(A 2 QCD /m c ) + 0(A QCD( ?/ml)) + 

1 /j, 2 q 2 

— -£■ — - + polynomial . (13) 

e 6 3 mj 

Here only the terms 0{(f) are kept. We also do not discuss perturbative corrections. 
Writing out the dispersion relation in e, 

lix J e — e 

\'hl rfg<4(e ~ f)+ r^I mw ^^ ft+h'hf dH2 < A & ^ +■■■ ^ 

and expanding it in 1/e we get the sum rules by equating the coefficients in front of 
l/e n . Here w AA = 2\mh AA . 
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3. Now let us discuss the phenomenological side of the sum rule. The structure 
function w AA is non-vanishing for positive e, 



12 



where the sum runs over all possible final hadronic states, the term with i — 
corresponds to the "elastic" transition B — > D* while i = 1, 2, . . . represent excited 
states with the energies E i = M i J r q 2 /(2M i ). Strictly speaking, |Fb_^| 2 does not 
present the square of a form factor; rather this is the contribution to the given 
structure function coming from the multiplet of the degenerate states which includes 
summation over spin states as well. In the particular example considered D is not 
produced in the elastic transition, so that in the elastic part one needs to sum 
only over polarization of D*. Therefore, the term "form factor" for Fs^i is rather 
symbolic. |Fb_>j| 2 depends on q. Moreover, Si in Eq. fll5| ) is the excitation energy 
(including the corresponding kinetic energy), 

Si = Ei — Ed* . 

For the elastic transition Sq vanishes, of course. 

The dispersion representation (|TJ]) and Eq. ( Jl3| ) lead to the following sum rule 



for the second moment of w AA (the coefficient in front of 1/e 3 , the third sum rule in 
the nomenclature of Ref . |]J ) : 

1 r °° I T? I 2 1 7& 

— dee 2 w AA {e) = V ^^^S 2 % = -/i 2 ^. (16) 

2vr J 1 W fr[ 2Ei l 3^m 2 c V ; 

We pause here to make a few remarks regarding Eq. (JIB]). First of all, since 
Sq = 0, this kills the elastic contribution in the left-hand side, and the sum actually 
starts from the first excitation. Second, since all Sf are of order Aq CD we need to 
know F B ^ 2 , Fb-*3, etc. only to the zero order in Aq CD . To this order all transition 
form factors to the excited states are proportional to q, i.e. 



'B- 



|2 „ -$, 



cx f. (17) 



(As a matter of fact, the transitions to P-wave states are relevant, see [Hj| for further 
details.) Moreover, taking account of Eq. ( |T7[ ) we can neglect 0(<f) part in S^s, so 
that in Eq. (|lj) 

Si = Mi-M D *. 

Third, m~ 2 in the right-hand side can be replaced, to the accuracy desired, by 
(M D *)~ 2 or by the mass of any excited state. 

After all these simplifications the third sum rule takes the form 

oo I rp 12 1 

g^(M,-M D .f = I^, (18) 



where 

Jt= 3- 

V M 

(it does not matter whose particular mass, Mb* or Mi, stands in the denominator). 
The next steps are rather obvious. The lower bound on /j% is a consequence of 
positivity of all individual contributions in the left-hand side of Eq. (0). Indeed, 
let us rewrite it as follows: 



1 oo ip 12 oo I rp_ 1 2 

3^ l ' f-i 2M i v 2 H 2M i v 2 



rf = «-gi^- + ESfS-W-<3- (i 9 ) 



The second term is evidently positive. The first sum can be found, in turn, by using 
the Bjorken sum rule ||16|| . This sum rule relates the sum over the P-wave states in 



the brackets to the (f dependence of the "elastic" B — > D* transition (the slope of 



the Isgur-Wise function |T7|). 



4. It is instructive to briefly reiterate derivation of the Bjorken sum rule, which, 
as explained above, is needed only in the zero order in Aqcd- Equating the coeffi- 
cients of 1/e in Eqs. ( |i~3"D and fli~4]) one immediately finds 

1 /cfc^(e) = £%?L = l * (20) 



2tt J 1 w £- 2£, : i 



The elastic part here can be parametrized in terms of the Isgur-Wise function S,{viv 2 ) 



17, 18 . The B — > D* transition has the form 



(D*{v 2 )\A^\B{ Vl )) = ^M B M D , [t^l + Vl v 2 ) - {ev 1 )v 2lx \i{v x v 2 ) 
where v 1,2 are the four-velocities. This means that 

(2£ D .r|FB- D .| 2 = ^ (i^) 2 If («.« 2 )l 2 - 1 - A' 2 ■ (2D 



Here p is the slope parameter [16] 



aviv 2 ) = 1 - p 2 ( Vl v 2 - 1) + ... = 1 - p 2 y + ... (22) 

and we used the fact that £ at zero recoil is unity M. 

Notice that although we discuss the Bjorken sum rule for the axial current actu- 
ally it can be derived for arbitrary current j a = bT a c. To the leading order in 1/m^c 
the universal factorization (|5|) takes place for the structure functions uw> = 21m h a i>. 
Moreover, the sum over any HQET degenerate multiplet of states gives 

' T,( B \ji\ H c)(Hc\j b \B) = C ab M Hc 1 -±^\U(viV 2 )\ 2 (23) 



2M, 



B 



where £# c is the Isgur-Wise function for the H c multiplet. 

At v = the sum rule ( p0[) is trivially satisfied since at zero recoil all inelastic 
form factors vanish, and we are left with the elastic contribution which reduces to 
unity. The term linear in v 2 yields a relation between the slope of £ and the inelastic 
contributions, 

1 ^ \F F 



P 



7 = E 



'B- 



(24) 



4 £J 2M i v 2 
Let us remind that the ratio \Fs^,i\ 2 /v 2 has the finite limit at zero recoil. Eq. (|24"D is 



the Bjorken sum rule proper fL6 |. Let us add for completeness that in the notations 
of Ref. fli~5"f , where the P wave inelastic contributions are written out explicitly, it 
takes the form 



p 2 -7 = EI-i ( /1(i)l 2 + 2£ 



lr 3 ( /1(l)l 2 



n=l n=l 

(for a simple derivation see Ref. ||13| ) . From these expressions it follows, in partic- 
ular, that p 2 > 1/4. 

Combining Eq. ( ]24j ) with Eq. ( |T9"D we finally arrive at 



V* 



36( p 



4 



+ 3E 



DO |p 

rB-n\ 



2M 4 v 



,72 



(tf-*?) 



(25) 



i=2 

8 i = M i - M D * . 



Eq. ( p5| ) is a direct n = 3 generalization of Voloshin's sum rule written for n = 2 
||, see also 0, 

1 _°°. I W I 2 

x = 2J ' l ' 2 -I» + «W l4 " i ' 1 ' (26) 

Since the second term in Eq. (|25| ) is positive we get the following obvious in- 
equality: 

/4 > 35 2 (p 2 - \) (27) 

(we remind that 5\ here is the lowest excitation energy, S\ = Mi — Mjj*). For a first, 
rough, estimate let us assume that S\ ~ 500 MeV and use the central value of the 
measured slope |B| of the B — > D* form factor for p 2 , 



Then 



0.84 ±0.12 ±0.08. 



pi > 0.45 GeV 2 . 



(28) 



Three comments are in order here regarding the sum rules presented above. 
First, the very same final results are obtained irrespectively of what currents we 
start from, axial or vector, or a mixture of these two. The only difference is that, 
say, for the vector currents we would get M D , not M D * in the definition of 5i. This 
difference is unimportant in the limit m\, )C — > oo, of course. This remark brings us to 



the second point. In Eq. fl2~7|) all subleading l/m bc terms have been omitted; thus, 
all quantities there refer to the infinite mass (static) limit, and the corresponding 
hadronic parameters must be understood just in this sense. In other words, rather 
than using the experimental value of p 2 and 5± measured in the beauty-to-charm 
transitions and the charmed family, respectively, one should use the static values of 
p 2 and 5i. (Hence, Mb is indeed equal to Mp. with our accuracy.) Finally, in the 
original sum rules the sum runs over all states including those which represent high- 
energy excitations described, in the dual sense, by perturbative formulae (see Ref. |I[ 
for more details). To get predictions for p 2 and A normalized at a low (quark- mass 
independent) scale p one must truncate the sum over the excited states at 5i ~ p 
and invoke duality between the perturbative corrections and the contributions of 
the excited states above p. 

In general, the sum rules at non-zero recoil get Aqcd/^&,c corrections which de- 
pend on the particular choice of the weak current considered and can be sizable. 
However, all corrections to the hadronic tensor h a b start with terms explicitly pro- 
portional to Aq CD /77t| )C 0, |2(J, ^TJ, see Eq. (A.l) in Ref. ||. The question is where 



the linear corrections come from? A source of subleading corrections is quite obvi- 
ous: they appear at the stage when one expresses e q in the theoretical formulae in 
terms of e; since Mb = m^ + A + ... (and the same for the charmed quark) they 
contain linear terms. This does not affect, of course, the first sum rule (n = 1), and 
in this case the prediction starts from unity plus corrections at the level A /m 2 c || . 

5. We proceed now to a more careful discussion of the numerical situation. The 
experimentally measured B —>■ D*(unpolarized) Iv decay rate is expressed in terms 
of the Isgur-Wise function in the leading approximation, see Eq. (|2l|) . In this 
approximation the slope of the Isgur-Wise function is related to the q 2 dependence 
of the B — > D* rate. It is clear that with l/m^^ and radiative corrections included 
the q 2 dependence of the decay rate does not exactly coincide any more with the 
slope of the Isgur-Wise function. These corrections were estimated in the literature 
(see the review paper f22[). Their effect seems to be equivalent to a decrease in the 
slope of the Isgur-Wise function, by about 9% . At the same time the radiative 
perturbative corrections were found [] to increase the slope by ~ 20%. Therefore, 
taking these estimates at their face value we are inclined to conclude that 

p 2 « p 2 exp - 0.1 (29) 

where p 2 in the left-hand side is the genuine static value of the slope parameter 
while p 2 is the value of the parameter obtained from experimental fits of the 
dependence of the B — > D*(unpolarized) Iv decay rate on the D* velocity. It is 
worth emphasizing again that the slope p 2 of the Isgur-Wise function by definition 
does not depend on the structure of weak currents considered. The above numerical 



1 Although the sign of the perturbative effects is obvious on general physical grounds, we think 
that the concrete procedure of evaluating them described in Ref. |22) systematically overestimates 
the velocity dependence, at least for the observable we discuss here. 
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estimates of both, perturbative and l/m& )C corrections, have been obtained for the 
real V — A current to which experimental numbers refer. 

Similar effects due to the finite mass of the c quark enter our lower bound im- 
plicitly when we use the observed mass values of the excited charmed mesons. In 
the future these pre-asymptotic corrections can be isolated in a model-independent 
way once the masses of the beauty counterparts are measured. The most sizable 
corrections are expected due to the chromomagnetic interaction of the heavy quark 
spin inducing hyperfine splitting among the members of the heavy spin multiplets. 
In particular, Mp* — Mp ~ 140 MeV. This effect is presumably accounted for by 
substituting the spin averaged masses for the ground S'-wave states and for the P- 
wave excitations, rather than actual masses of D, D*, etc. We actually did this spin 
averaging. Another shift arises due to the heavy quark kinetic energy term in the 
hadron mass. It is natural to expect its value to be smaller in the excited mesons 
than for the ground state. Therefore, the static limit of 5% is expected to be some- 
what larger than the value of S\ experimentally observed for the actual charmed 
particles, but probably not more than by 50 MeV. We then use the value 

<Ji « 500 MeV (30) 

as a very reasonable educated guess for the static value of S\. 
With the parameters from Eqs. (EST), (BUD we fhiahy get 



& > 0.37 GeV 2 , 

A > 500 MeV (31) 

where the second relation comes from Voloshin's sum rule (|26|). These lower bounds 
are seen to lie not very far from the estimates obtained earlier within QCD sum 



rules 23 



(4 ~ 0.55 GeV 2 , A ~ 450 MeV . (32) 

Note that the lower bound on /z 2 in Eq. (RTI) is numerically quite close to the bound 



(HI) derived recently in || |!|. 

Unfortunately, numerical uncertainties in all the numbers above prevent us from 
making too strong a statement. Nevertheless, let us assume for a moment that 
future refined measurements and calculations of the subleading corrections in the 
third sum rule will confirm these values and establish the fact that two inequalities 
in Eq. ( |3lD are rather close to saturation. This would mean that the sum rules are 
actually saturated - to a reasonable degree of accuracy - by the contributions from 
the states with masses around M D + Si generically called D** in this context. To 
account for nonperturbative effects in b — > c decays one then would need only to 
consider one inelastic channel, "_D**" . The higher excited states will be represented 
(in the sense of duality) by purely perturbative probabilities calculated in the free 
quark-gluon approximation. We actually consider such a situation as a most natural 
scenario in QCD. It is worth noting that the Dtt contribution to the third sum rule 

9 



is suppressed for soft pions, unlike the first sum rule, where it was quite substantial 
0. The effect of the "hard" pion emission is well represented by some of the P-wave 
D** resonances. 

6. We have derived the third sum rule for the b — > c transition in the SV limit 
and showed how one can use it to constrain the kinetic energy parameter /j% by using 
the data on B — ► D* . In principle it is quite conceivable that the full differential 
distribution in q and (f in the inclusive semileptonic B decays will be measured in 
the future. This measurement can then be immediately translated in the value of 
fj%i one °f t ne most important parameters of the heavy quark physics. The more one 
will learn about the decays to the excited states the more accurate the determination 
of \x^ will become. 
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